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ABSTRACT 

Aims. IKT 16 is an X-ray and radio-faint supernova remnant (SNR) in the Small Magellanic Cloud (SMC). A detailed X-ray study 
of this SNR with XMM-Newton confirmed the presence of a hard X-ray source near its centre, indicating the detection of the first 
composite SNR in the SMC. With a dedicated Chandra observation we aim to resolve the point source and confirm its nature. We 
also acquire new ATCA observations of the source at 2.1 GHz with improved flux density estimates and resolution. 

Methods. We perform detailed spatial and spectral analysis of the source. With the highest resolution X-ray and radio image of the 
centre of the SNR available today, we resolve the source and confirm its pulsar wind nebula (PWN) nature. Further, we constrain the 
geometrical parameters of the PWN and perform spectral analysis for the point source and the PWN separately. We also test for the 
radial variations of the PWN spectrum and its possible east west asymmetry. 

Results. The X-ray source at the centre of IKT 16 can be resolved into a symmetrical elongated feature centring a point source, the 
putative pulsar. Spatial modelling indicates an extent of 5.2"of the feature with its axis inclined at 82° east from north, aligned with a 
larger radio feature consisting of two lobes almost symmetrical about the X-ray source. The picture is consistent with a PWN which 
has not yet collided with the reverse shock. The point source is about three times brighter than the PWN and has a hard spectrum of 
spectral index 1.1 compared to a value 2.2 for the PWN. This points to the presence of a pulsar dominated by non-thermal emission. 
The expected E is ~ 10 37 erg s~' and spin period ; 100 ms. However, the presence of a compact nebula unresolved by Chandra at the 
distance of the SMC cannot completely be ruled out. 

Key words, galaxies: Magellanic Clouds - ISM: supernova remnants 


1. Introduction 

Composite supernova remnants (SNRs) are a subclass of core 
collapse SNRs where non-thermal radio emission is observed 
from both the expanding shell of the SNR, and from the PWN 
located inside it. The PWNe are powered by the relativistic out¬ 
flows from the young neutron star and can be observed at all 
wavelengths, but mostly in the radio and X-rays. The PWN mor¬ 
phology can provide crucial information on the properties of the 
outflow, the interacting a mbient medium and the geometry of the 
pulsar powering it (see lGaensler & Sland l2006. and references 
therein for detailed understanding of the structure and evolution 
of PWNe). The Chandra observatory with its excellent spatial 
resolution and high sensitivity has provided significant break¬ 
throughs in the study of PWNe. Apart from steadily increasing 
the census on the number of PWNe discovered, it has provided 
a unique opportunity to study in detail their spatial and spec¬ 
tral struct ures and signatures of int eraction with the surrounding 
medium dKargaltsev & Pavlovi2008l) . Most of these are however 
Galactic sources and there are no confirmed PWNe in our satel¬ 
lite galaxy the SMC, yet. To this day a_ total of 24 classified 
SNRs are known in the SMC (lHaberl et al .1 120 1 2bl and refer¬ 
ences therein). Out of them very few are PWN candidates. One 
such source, HFPK 334 has_been recently dismissed as a back¬ 
ground source (ICrawford et al .1 [2014l) . Discovery of a PWN in 
the SMC will open a new window in the study of rotation pow¬ 


ered pulsars in the satellite galaxy which has a rich history of 
active star formation and is therefore expected to host young en¬ 
ergetic pulsars. 

IKT 16 is an X-ray and radio-faint SNR in the SMC first 
studie d with XMM-Newton in a su rvey of known SNRs in the 
SMC (Ivan der Hevden et alJ 120041) . In this study a region of 
harder X-ray emission was found at the centre of the remnant al¬ 
though it could not be probed further due to poor statistics with 
a single exposure. [Owen et akl d201 lb further used eight more 
archival XMM-Newton observations (total useful exposure 125 
ks and off-axis angle 8-12') taken subsequently to carry out a 
multiwavelength study of the spatial and spectral properties of 
this SNR and the associated central source. The authors found 
substantial evidence that the unresolved source detected at the 
centre of the SNR by XMM-Newton, is a PWN associated with 
it. Radio images from the ATCA and MOST surveys displayed 
faint radio structures correlated with the X-ray throughout the 
remnant. The brightest feature was an extended radio emission 
corresponding to the X-ray source and extending a distance of 
40" towards the centre of the SNR dOwen et alJ 1201 it) . This 
picture is consistent with a moving pulsar in the SNR. Further, 
the large size of the remnant suggested that the SNR is in the 
adiabatic Sedov phase of evolution. In the Sedov model the age 
of the SNR was ~ 14.7 kyr and implied the PWN may have inter¬ 
acted with and been compressed by the reverse shock. The bright 
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central source was located 30" from the SNR centre implying a 
transverse kick velocity of ~ 580 km s . 

To further resolve the nature of the central source and es¬ 
tablish the presence of the first PWN in the SMC, a 40 ks on- 
axis Chandra observation was solicited. In this paper we de¬ 
scribe the results of the Chandra observation of IKT 16. We have 
performed detailed spatial and spectral analysis of the hard X- 
ray emission near the centre of the SNR, previously unresolved. 
With the unprecedented spatial resolution of Chandra we have 
resolved the source into a central point source (a putative pulsar) 
and an extended emission surrounding it, indicating a PWN na¬ 
ture for the source. The observations and analysis are described 
in section 2. Section 3 presents the spatial analysis including 
the imaging and morphological fitting of the source. Section 
4 presents the detailed spectral analysis of the central source, 
the surrounding nebula and its decomposed regions. Section 5 
presents the discussion and section 6 the conclusions. 


2. Observations & Analysis 

2.1. Chandra ACIS observation and data reduction 

The Chandra observation (ObsID 13773) was carried out with 
the ACIS-S as the primary instrument in the timed exposure 
mode. The observation was performed on 09-02-2013 with the 
central source of IKT 16 positioned at the aimpoint of the S3 
CCD with an exposure of 38.5 ks. It is worthwhile to mention 
here that there was a previous Chandra observation (obsID 2948, 
9 ks) where the source was detected (lEvans et al.ll2010l) . But the 
short duration of this observation and its large off axis angle pre¬ 
vented further investigation of the source from this data. Data 
reduction was performed using CIAO 4.6 using CALDB 4.6.3 
and the standard analysis procedure prescribed Q. The level 1 
event file was reprocessed with chandraj-epro , which incorpo¬ 
rates the subpixel repositioning algorithm EDSER as a default 
for attaining better angular resolution for sources near the centre 
of the FOV. The effective exposure of the observation after filter¬ 
ing was 38.5 ks. We also checked for possible presence of pileup 
in the data. Using the tool pileup jnap, the estimated fraction of 
pileup in the centremost pixel is ; 5%. Hence the effect is not 
important for our observation . 

Images were created using the task dmcopy. Spatial analy¬ 
sis was performed using the Sherpa analysis package 4.40. The 
task specextract was used for extracting source and background 
spectra and response files from regions of interest. The selection 
of regions used for spectral extraction is described in the spectral 
analysis section. Spectra were fitted using XSPEC vl2.8. 

2.2. Imaging 

Figure [Qi shows a full resolution Chandra image (0.5-8 keV) 
centred on the source near the centre of IKT 16, and zoomed into 
a region 20"x 25" wide. The source appears to be symmetrically 
elongated in the east-west direction with a bright point source 
located at its centre. The elongation measures about 5". To ac¬ 
curately determine the position of the central source, we created 
a subpixel image of the same (at ~ of the ACIS pixel resolu¬ 
tion) and applied the source detection algorithm celldetect. The 
coordinates of the point source are RA(J2000)=00 / '58'"16.85 i 2 3 


1 http://cxc.harvard.edu/ciao/ 

2 http://cxc.harvai'd.edu/sherpa4.4/ 

3 http://heasarc.gsfc.nasa.gov/xanadu/xspec/ 


Dec=-72°18' 05.60" considering an error of 0.6" at 90 % con¬ 
fidence level in absolute Chandra astrometry. The net count rate 
from the point source after subtracting the nebular component is 
0.011 c/s, and from the entire nebula after subtracting the point 
source contribution is 0.003 c/s. The details of background sub¬ 
traction are discussed in the spectral analysis section[4] The point 
source is about three times brighter than the extended emission. 
After removing the contribution from the bright point source at 
the centre, and smoothing the image with a Gaussian kernel of 
width 5", there is evidence of a diffuse emission extending fur¬ 
ther out up to ~ 30" (see Fig. |TJl). This diffuse component is 
discussed in Sect l4.5l The same image shows hints of the pres¬ 
ence of the SNR, especially the excess in the north coincident 
with that in the XMM - Newton image of IKT 16 (i Owen et a li 

l20ll . 


2.3. Radio observations 

Our new ATCA observations, project C2521 (Cl: J. van Loon), 
used the Compact Array Broadband Backend (CABB) with the 
6A array configuration at 2.1 GHz providing improved flux den¬ 
sity estimates and resolution. These new images were acquired 
on 02-01-2012 with ~ 10.48 hours integration over the 12 hour 
observing session. The radio galaxy PKS 1934-638 was used as 
a primary flux calibration source for all observations, with the 
radio sources PKS 0230-790 and PKS 2353-686 used for phase 
calibration. A standard calibration process was carried out using 
the miriad data reduction software package (Sault et al. 2011). 
In order to improve the fidelity and sensitivity of the final im¬ 
age, a single iteration of self-calibration was performed on the 
strongest sources in the field. A uniform weighting scheme was 
subsequently used throughout the imaging process, as it pro¬ 
vided a balance between improving the theoretical rms noise 
while maintaining an adequate beam shape. Given the 2 GHz 
of bandwidth provided by CABB, images were formed using 
miriad multi-frequency synthesis (mfclean; Sault & Wieringa 
1994). The same procedure was used for both U and Q Stokes 
parameter maps. However, there was no reliable detection in the 
U or Q intensity parameters associated with this object, imply¬ 
ing a lack of polarisation at lower radio frequencies (2.1 GHz). 
The final image produced (after primary beam corrected), which 
possess a FWHM of 3.38" x 3.16" and a PA of -17.8°, has a lcr 
rms noise level of 16ply. In Fig.[2]we show the ATCA 2.1 GHz 
surface brightness image of PWN IKT 16. 


3. X-ray Spatial analysis 

Taking advantage of the excellent spatial resolution of Chandra, 
we performed detailed spatial analysis of IKT 16 for the first 
time. To provide further evidence of the extended emission, and 
to probe the source extent we created a radial profile centred on 
the point source from the data. We compared this with the two 
dimensional model of the point spread function (PSF) to look for 
an excess indicating the extended emission. Further, we tested 
for signatures of asymmetry along the two halves of the nebular 
emission. Finally we performed a morphological fitting of the 
elongated structure with a simple model of the nebula centred 
on the point source, and determined its geometrical parameters 
like the size, ellipticity and the rotation angle. 
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Fig. 1 . a (top left)'. Chandra ACIS-S (0.5-8 keV) full resolution image of the PWN in IKT 16. The image size is 20"x 25". The scale 
is in square root, the X and Y axis in degrees, and units are counts for all the images. The brightest pixel in the centre corresponds 
to the putative pulsar, b (top right): Best-fit model of the PWN (0.5-8 keV). Overlayed are contours from the data which have been 
smoothed with a Gaussian of cr = 3". The contour levels are plotted at values of 0.7, 1.6 and 30 counts arcsecs 2 . c (bottom left): 
Same as the top left figure showing the regions used for spectral extraction. The central circle corresponds to the pulsar, and the 
rectangular box the entire nebula with the central point source removed. The dashed circle corresponds to the outer boundary of the 
inner, and the inner boundary of the outer nebular extraction regions respectively, d (bottom right): Larger (0.5-2 keV) image centred 
on the PWN of IKT 16 (ACIS-S3), with the bottom right corner of the image corresponding to ACIS-S2. The box region used for 
nebular extraction is shown in black solid lines and the background annular region used for spectral extraction (inner radius 34") is 
shown in white solid lines. The white dashed circle indicates the position and extent of the SNR in lOwen et alJ (1201 ll) . The point 
source contribution at the centre has been removed and replaced with values obtained by interpolation from the nebular region. The 
image has been smoothed using a Gaussian kernel of width 5". 

3.1. Radial profile 


Radial profiles were created up to 1' in two energy bands (0.5- 
2 keV, and 2-7 keV) by extracting net counts in circular annuli 
centred on the point source using the tool dmextract and then 
rebinning to reach a reasonable statistical precision. The back¬ 
ground was extracted from a circular region of the same area, 
away from the source region but inside the SNR. The PSF of the 
observation was simulated using the Chandra ray tracer chaRT 


0 which simulates the High Resolution Mirror Assembly based 
on the energy spectrum of the source and the observation ex¬ 
posure. The output of chaRT was modelled with the software 
MARX 0 taking into account the instrumental effects and the 
EDSER subpixel algorithm to be consistent with the observa¬ 
tional data. The best-fit spectrum of the point source (see section 
ED was used. Fig. [3] shows the radial profile of the observa¬ 
tion along with the simulated PSF. The data at both the energy 


4 http://cxc.harvard.edu/chart/ 

5 http://cxc.harvard.edu/chart/threads/marx/ 
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Fig. 2. ATCA 2.1 GHz image of PWN IKT 16 overlaid with Chandra contours (1.25, 5, 25 and 100 c/arcsec 2 ). The cross marked in 
green denotes the centre of the SNR in IKT 16. The synthesis beam size of 3.38" x 3.16" at PA of -17.8° is shown on the bottom 
left corner. The 2.1 GHz image grayscale bar is shown in units of mJy/beam. 


bands are consistent with a point source up to a radius j 1" be¬ 
yond which it clearly has higher net counts than expected from 
a point source simulated at the same position and with the same 
spectral parameters. This corresponds to the nebular component. 
Although the geometrical model (see section [131) indicates that 
the source FWHM is about 5", this exercise indicates that the 
source extends further beyond. This excess is best seen in the en¬ 
ergy range of 0.5-2 keV. Its brightness decreases outwards like 
r L5 . Figure [lji indicates that this further extension is probably 
elongated in the east-west direction like the main X-ray nebula 
and the radio nebula. 

3.2. Investigating signatures of asymmetry in the nebula 
morphology 

Visual inspection of the extended source near the centre of IKT 
16 from Fig. |T^, gives an impression of an elongated symmetric 
structure. However, we looked for the possibility of an east-west 
asymmetry by dividing the nebular structure into two halves 
along the axis going through the point source, and compared 
the total counts from the two regions. The number of counts in 
the east and west nebula are 50 + 9 and 67+10 respectively 
indicating that the west nebula may be slightly brighter than its 
other half. It cannot be determined conclusively as the values are 
consistent within their errors. We investigated this issue further 
through spectral analysis of the same regions in section l4~4l 

3.3. Spatial modelling 

Spatial fitting was performed on a larger region of size 35"x 
45"centring the source to detect the existence of a possible dif¬ 
fuse emission underlying the nebular one, and also constrain the 
background better. The PSF described in the previous section 
was loaded as a table model in Sherpa to model the point source 


Table 1. Parameters of the best-fit model to pulsar wind neb¬ 
ula in IKT 16. Model used is a 2-D Gaussian function for the 
nebula and a constant background. Theta is the angle between 
the major axis and the north direction. 


Parameter 

Value 

Units 

FWHM 

c 9-u.y 

J +10 

arcsec 

Ellipticity 

0.6 ±0.1 


Theta 

82 ±7 

degrees 

Amplitude 

2.9 ±1 

counts/pixel 

Background 

0.12 ±0.01 

counts/pixel 


emission. The remaining excess consisted of the background and 
the emission from the nebula. This was modelled with a constant 
background and a 2-D unnormalized Gaussian function (gauss2d 
model in Sherpa), convolved with the PSF. In the initial itera¬ 
tions it was found that the position of the 2-D Gaussian used 
to model the nebula was consistent with the point source. In or¬ 
der to constrain the geometrical parameters of the nebula like its 
FWHM along its major axis, ellipticity and rotation angle better, 
its position was fixed to the point source henceforth. The best-fit 
parameters were determined by the C-statistic (iCas I) 1979 ) and 
errors were estimated at 90% confidence level. The residuals do 
not show any systematic pattern indicating that this analysis is 
not sensitive to any substructures apart from the elongated struc¬ 
ture of the nebula. Fig. [TJt shows the model with the contours 
from the data smoothed with a Gaussian of cr — 3"overlayed 
on it. This highlights the fact that the model is a good descrip¬ 
tion of the data. It is worth mentioning that we also tried to fit 
a 2-D Lorentz model with a varying power-law (beta2d model 
in Sherpa) instead of the 2-D Gaussian used to model the neb¬ 
ula. We noticed that by freezing the rotation angle and ellipticity, 
the power-law index a tends towards the maximum limit of 10 
which closely approximates the 2-D Gaussian model. The dif- 
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Radius [arcseconds] 

Fig. 3. Radial counts profile of the source (blue star: 0.5-2 keV; 
red triangle: 2-7 keV) plotted against the 2-D PSF (0.5-8 keV, 
in black crosses), clearly showing the extended nature of the 
source. 


ference in C-statistic between beta2d and gauss2d is 15 in the 
Lorentzian limit (a -1 in beta2d), and 2 in the Gaussian limit 
(a=10 in beta2d). We therefore concluded that the data does not 
favor broad wings and a 2-D Gaussian is preferred for the nebu¬ 
lar emission. 

The best fit parameters are listed in Table Q] The estimated 
FWHM of 5.2" is in agreement with that measured from the im¬ 
age of the source. It is inclined at an angle of 82° between the 
long axis and the north direction and is aligned with the radio 
nebula (see Fig. [2}. An ellipticity of 0.6 indicates a major to mi¬ 
nor axis ratio of 0.8. Another important result is that there is 
no signature of displacement between the point source and the 
nebula. The best-fit position of the centre of the nebula obtained 
corresponds to R.A.(J2000)=00:58:16.824 and DEC.(J2000)=- 
72.18:05.32 with an error of 0.1"and 0.2"on the x and y po¬ 
sitions respectively. This is consistent with the position of the 
point source obtained with celldetect. 


4. X-ray Spectral analysis 

We have performed detailed spectral analysis of the point source 
and the nebular component of IKT 16. In addition we have also 
investigated the diffuse component extending to larger scales. 
The main deciding factors for the spectral analysis are the re¬ 
gions used for the spectral extraction and the background mod¬ 
elling, taking into account the contribution of the point source 
in the nebular spectrum and vice versa. This is described in the 
subsequent subsections. The analysis was performed in the en¬ 
ergy range of 0.5-7 keV. C-statistic was used for spectral fit¬ 
ting and errors were estimated at 90% confidence interval. To 
account for the photoelectric absorption by the interstellar gas, 
two absorption components were used as in lOwen et al.l (1201 ll) . 
The first, phab s component was fixed at the Galactic value of 
6xl0 20 erg s _1 (Dickey & Lockmanl 199(j) . and a free absorption 
(v phabs) component was used to account for absorption inside 
the SMC. This second component has metal abundances fix ed at 
0.2 solar, as is typical in the SMC iRussell & Dopitall 19921) . 

4. 1. Point source spectrum 

For the point source, a circular region of radius 1.5 pixels was 
extracted centred on the best-fit coordinates of the source. The 
extraction region is shown in Fig.Q]:. A circular annulus outside 


the nebular extraction region was used for the background spec¬ 
trum. Apart from this the astrophysical background due to the 
nebular contribution was also considered. For this we calculated 
the fraction of nebular emission contributing in the point source 
extraction region (from the best-fit morphological model of the 
source), and accounted for it as an additional model component 
in the spectral fitting, with parameters fixed to the best fit values 
obtained in section l4~2l 

We tried to fit the spectrum with several models including an 
absorbed power-law, an absorbed blackbody and a combination 
of both models. The absorbed blackbody model leads to an un¬ 
acceptable spectral fit with a difference of 54 in C-statistic value 
between the two models with respect to the power-law model. 
The absorbed power-law model provides a good fit as illustrated 
on Fig. [4jt. The best fit parameters along with the flux are tab¬ 
ulated in Table [2] The addition of a blackbody component was 
not required and did not improve the fit significantly. Also, the 
additional absorption column density inside the SMC could not 
be constrained well and was prone to large error bars as can be 
seen from Table[2] Its value was however consistent with that ob¬ 
tained from the previous analysis of the XMM-Newton observa¬ 
tions bv lOwen et alJ ( 2011 1 which had smaller errors associated 
with it probably due to the larger number of counts. 

4.2. Nebula spectrum 

For the spectral extraction of the extended source, a box region 
was used with the region corresponding to the point source ex¬ 
cised (see Fig. |TJ:). Its size was optimized from the morphologi¬ 
cal modelling of the nebula, and its centre was made coincident 
with the point source. The background was taken from the same 
annular region that was used as a background for the point source 
spectral extraction. For the additional astrophysical background 
contributed by the point source leakage in the nebular region, 
we extracted a spectrum from the events from the PSF simulated 
by chaRT with spectral parameters similar to those of the point 
source, in the same region that was used for the nebular spectral 
extraction, and added it to the internal background spectrum as 
an additional background component. 

We used an absorbed power-law as the spectral model for 
all the nebular fits described henceforth. In contrast to the 
point source spectrum, the spectrum of the nebula is very soft. 
Keeping the absorption inside the SMC free in these fits leads 
to strong correlation between the absorption column density and 
the power-law index T. To avoid this it is better to constrain this 
parameter from the point source spectrum, and freeze it to this 
value for the nebular spectral fits. However, the better capability 
of the XMM observations of IKT 16 to constrain the local ab¬ 
sorption density, as discussed earlier, led us to fix this value to 
that obtained from the previous analysis of IKT 16 with XMM- 
Newton dOwen et alJl201 lh in all the nebular fits. The spectrum 
along with its best-fit model and residuals is shown in Fig. [4j>, 
and the best fit parameters are tabulated in Table [2] 

4.3. Outer and Inner nebula spectra 

In order to look for changes in the spectral parameters, particu¬ 
larly a spectral steepening with radius in the nebula, we divided 
the total nebular region into two parts, by choosing an inner an¬ 
nular region centred on the point source, and an outside region 
excluding it. The extraction regions are shown in Fig. Ut. The 
Fig. [5ji and [fiji show the inner and outer spectrum respectively 
with their best-fit models and residuals, and Table [2] their best- 
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Energy (keV) Energy (keV) 

(a) (b) 

Fig. 4. The figures show the ACIS-S spectra for the putative pulsar (a ,left) and the entire PWN (b ,righf). The solid lines correspond 
to the respective best-fit spectral models. The plots have been rebinned for visual clarity. 


fit parameters. Although the outer nebula is hinted to be softer 
then its inner counterpart by the absence of counts above 5 keV, 
we detect no steepening of the spectral index T which would be 
indicative of synchrotron cooling. The T values are consistent 
within errors at 90 % confidence. 

4.4. East-west nebula spectra 

In Sect. 13.21 where we report the extracted total number of counts 
in the east and west nebular region, we find evidence that the 
west of the nebula may be slightly brighter although the val¬ 
ues are consistent within errors. In continuation, we also extract 
spectra from the same regions following the same procedure as 
in Secs. 14.21 and 14.31 Figures [5}; and Fig. [5jl show the east and 
west spectrum respectively with their best-fit models and resid¬ 
uals, and Table [2] their best-fit parameters. As indicated by the 
earlier exercise in Sect. 13.21 the spectrum from the west of the 
nebula shows slightly higher value of flux but it is comparable to 
that measured from the east nebula within errors at 90 % confi¬ 
dence level. We do not detect any significant change in the spec¬ 
tral index F. 

4.5. Diffuse emission 

From Figs, [ljl and [3] we see evidence that the source extends 
further beyond the region adopted to study the nebula spectrum. 
To investigate further, we extracted the spectrum of the diffuse 
component by choosing an annular region excluding the nebu¬ 
lar extraction region (box region in Fig. [T]:) and extending up 
to the size of the XMM-Newton point source extraction region 
of 20" dQwen et al.ll201 ll) . Using the same spectral model as that 
for the nebular spectrum, we find that the diffuse component has 
a flux and spectral index (tabulated in Table 0 comparable to 
that of the main nebular component. We also tried to fit the spec¬ 
trum with a Sedov model, since that emission might be associ¬ 
ated with the SNR. We fixed the parameters to those obtained 
by lOwen et aD ( 201 1!) (the data quality does not allow fitting 
anything else than the normalisation). The resulting C-statistic 
is higher by 9 than that of the power-law model. This does not 
favor the thermal model, but does not allow ruling it out either. 
However, no other part of the SNR is as bright, and this diffuse 
emission is centred on the main X-ray nebula. So we favor the 
interpretation as a further PWN component. 


In addition, we also investigated whether this emission could 
be a halo due to the scattering from the foreg round dust in 
the interstellar medium. iPredehl & Schmitt] (11995b using ROSAT 
observations, studied X-ray scattering halos around 25 point 
sources including the three bright sources in the LMC, LMC X- 
1, LMC X-2 and LMC X-3. The local absorbing density around 
IKT 16 is lower that in LMC X-l, but higher than in LMC X-2 
and LMC X-3, and thus lies in the range covered in this study. 
The authors found that the relative intensity of the scattering ha¬ 
los of the LMC sources w.r.t the point source are ~ 1%. This is 
much weaker than what we find for the diffuse component which 
is ~ 25% (j 2 keV) of the putative pulsar emission. Moreover, 
ROSAT observed a flat profile for the halos extending to i 100". 
Within our extraction region of the diffuse emission (20"), the 
fraction is expected to be even lower than 1%. Hence the dust 
scattering origin for this diffuse emission is very unlikely. 

Finally, we compared our results with the spectral model of 
the point source seen with XMM-Newton (Table 3 in lOwen et al.l 
1201 ll . counting the spillover of the point source into the SNR 
extraction region, and the SNR emission in the point source re¬ 
gion). The total flux measured in the energy range of 0.5-8 keV 
from the Chandra observation is 21.3 x 10~ 14 erg cm~ 2 s _1 (pu¬ 
tative pulsar + nebula + diffuse emission; see Table [2] for the 
obtained values) comparable to 24.9 x 10 14 erg cm -2 s _1 mea¬ 
sured from the XMM-Newton spectral model. The unabsorbed 
luminosity in the same energy range (0.5-8 keV) measured from 
Chandra is 1.0 x 10 35 erg s 1 . In comparison, the XMM spectral 
model corresponds to an unabsorbed luminosity of 1.1 and 1.3 
xlO 35 erg s~ 1 in the 0.5-8 keV and 0.5-10 keV energy ranges 
and is in agreement with the Chandra results. It should be noted 
that, although the total unabsorbed luminosity (0.5-10 keV, SNR 
+ point source) in lOwen et al.l (1201 ll) is consistent with the spec¬ 
tral model, there seems to be an error in quoting the unabsorbed 
luminosity of the point source as 1.6 xlO 35 erg s _1 in that paper. 

5. Radio 

The PWN in the radio-continuum regime extends far beyond X- 
ray detection as can be seen in Fig. [2] The radio extent of the 
source as measured visually, is 70 + 10". We note that the X- 
ray point source correlates to one of the peaks in our RC image. 
However, no corresponding point source can be drawn anywhere 
near the X-ray point source. 


6 
















































C. Maitra et al.: IKT 16: the first X-ray confirmed composite SNR in the SMC 



Fig. 5. Same as in Fig. |4]but for the inner nebula a, outer nebula b, east nebula c and west nebula d. Please note that in the outer 
nebula, i.e. Fig. b, there are no events above 5 keV. 


Table 2. Parameters of the best-fit spectral model to putative pulsar and the nebula. Errors are quoted at 90 % confidence. The 
Chandra energy range is 0.5-8 keV. 


Region 

SMC absorption 
10 21 cm -2 

Power-law photon Index 

r 

Fx" 

0.5-8 keV 

Fx" 

2-10 keV 

Lx 

0.5-8 keV 

Putative pulsar 
Nebula 

Inner nebula 

Outer nebula 

West nebula 

East nebula 
Diffuse emission 

-3.5 

3.4 (f) c 

3.4 (f) c 

3.4 (f) c 

3.4 (ff 

3.4 (f) c 

3.4 (f) c 

1.11 ±0.23 

9 91 +0.40 
^•^-0.37 

9 99+0.48 

y -Q.4Q 

9 90+O.67 
-0.47 

2 ^+ 038 

_ Z - U -0.30_ 

16+1 

2.2±0.2 

1.6+0.3 

0 75 +038 

1.4+0.3 

0 93 +0 - 4 

u- 7 -0.36 

3.1+0.2 

17+2 

1 9+0.4 

10^ 

1- -0.7 

0.47±0.63 

0 9 +05 
0.55+0.67 

1 0+0.3 
-0.2 

7.20 

1.19 

0.78 

0.41 

0.67 

0.45 

1.70 


“ - Observed flux in units of 10~ 14 erg cur 2 s -1 
b - absorption corrected luminosity in units of 10 j4 erg s -1 assuming a distance of 61 kpc. 

(f) - Parameter fixed for consistency between fit regions. 

c - SMC absorption for all nebular spectra was frozen to the column density obtained from the XMM-Newton observation . 


We point to the new high resolution data (Fig.[2]i where two 
lobe like features, somewhat symmetric, each 20 + 5" in extent, 
outside of the X-ray nebula are seen. These radio-continuum fea¬ 
tures oppose our initial scenario that the radio nebula has been 
pushed aside by the reverse shock dOwen et al.ll201 il l. It seems 
it is more likely that we are still seeing the nebular expansion in 
the cold ejecta. 

Across the 2.1 GHz band, the radio spectral index is flat as 
reported in iOwen et al. (20JJ ). Somewhat surprising, we did not 


detect any polarisation in either Q, U or V stokes parameters 
even that our detection level is better than 1%. The small scale 
structure may play an important role in smoothing out any weak 
polarisation. At the same time, depolarisation may reveal re¬ 
gions of underlying turbulence and/or compression and heating 
of th ermal material at va rious shocks within the remnant sys¬ 
tem (lAnderson et al.l[l995l) . This depolarisation may also indi¬ 
cate somewhat older age of the remnant. The noticeably younger 
PWNe in the LMC such as N157B (lLazendic et al.l lToOO ). SNR 
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J0453 -6829 (lHaberl et alJ2012ah . and 0540-69.3 dBrantseg et al.1 
|20 1 4l) all have strong polarisation. 


6. Discussion 

Confirmation of the PWN nature of the source makes IKT 16 the 
first composite SNR to be discovered in the SMC. The nebular 
structure is consistent with the presence of an elongated PWN 
around a putative pulsar, which extends further out in the form 
of a fainter diffuse emission. However, we cannot rule out the 
presence of an inner compact X-ray nebula of the order of 0.1 pc 
which cannot be resolved at the distance of the SMC. Therefore, 
the overall impression is that of the bright central point source 
corresponding to a putative pulsar or a putative pulsar + unre¬ 
solved inner PWN component, and an extended emission com¬ 
ponent representing the PWN. 

Although the present Chandra observation in full-frame 
mode is not capable of detecting pulsations from the central pu¬ 
tative pulsar, we can nevertheless infer its properties to a cer¬ 
tain extent by assuming some basic characteristics of the neutron 
star as is described below. In the later subsections, we also dis¬ 
cuss the evolutionary stage of the pulsar wind nebula as expected 
from its morphology. 

6.1. Derived properties of the pulsar 

The nebular emission in X-rays reflects the youngest generation 
of the emitting particles. Therefore, a strong correlation is ex¬ 
pected between the non-thermal X-ray luminosity of the PWN 
and its pulsar, with the spin down properties of the pulsar it¬ 
self. The properties of the putati ve pulsar in IKT 16 were es¬ 
timated adopting the work of iKargaltsev & Pavlov! (1 20081 ). In 
this paper the correlation between the luminosity of the PWN, 
the non-thermal luminosity of the pulsar, and its various proper¬ 
ties are derived from a large sample of Chandra observations of 
PWNe. The unabsorbed luminosities of the PWN and the pul¬ 
sar in IKT 16 are reported in Table [2] Using these values and 
the tab ulated properties of the PWNe from Kargaltsev & Pavlovl 
d2008h . we estimated the spin-down power of the putative pul¬ 
sar by comparing the total non-thermal luminosity of the sources 
(PWN+pulsar component) with their spin-down power. This en¬ 
sured our estimate was free from the assumption that the central 
source is exclusively the pulsar and does not contain a compact 
nebular component. It is noteworthy to mention that in order to 
perform this estimation, we have not accounted for the flux from 
the diffuse emission (described in section 1431 ) which is not as 
well measured and whose origin is less unambiguous. At the 
same time addition of this component would not alter our results 
in the limits of the correlation measured in Karg altsev & Pavlovl 
(2008). Fig. [ 6 ] plots the spin-dow n power against the total lu - 
minosity for all the objects from IKargaltsev & Pavlovl ((2008), 
with a line drawn to indicate the total luminosity of the cen¬ 
tral source in IKT 16. It corresponds to an expected spin-down 
power E ~ 10 17 erg s' 1 . Further, a ssuming its age to be the same 
as that measured from the SNR (lOwen et al .112011l) . i.e. 14.7 
kyr, we overlayed th e central source on the P - P diagram from 
the ATNF catalogue dManchester et al.ll2005h . The expected spin 
period is j 100 ms with P ~ 10' 13 s s' 1 . This points towards a 
young and energetic pulsar powering the PWN in IKT 16, in fact 
the youngest in the SMC. Future on-axis XMM observations or 
Chandra observations using operating modes with better timing 
resolution can be used to search for pulsations from this source. 



Log E aol (0.5-8 keV) 

Fig. 6. Dependence of the total non-thermal luminosity (neb- 
ula+pulsar) on the pulsar spin-down power for the PWNe 
observ ed with Chandra . Adapted from Kar galtsev & Pavlovl 
( 2008 ). The red dashed line indicates the luminosity of the puta¬ 
tive pulsar + PWN in IKT 16. 


6.2. Energy spectrum 

The spectra of both the central and the extended source can be 
satisfactorily fitted with a power-law model, expected from a 
synchrotron emission dominated spectrum. The hard power-law 
spectral index 1.1 for the point source resembles that of a typical 
young rotation powered pulsar where the non-thermal radiation 
is generated by the particles accelerated in the pulsar magneto¬ 
sphere dCurtisll992l and references therein). No additional ther¬ 
mal component (which would denote the emission from the sur¬ 
face or the polar caps of the neutron star) was required to fit the 
pulsar spectrum. The point source is also three times brighter 
than the surrounding nebula. If this comprises exclusively the 
pulsar emission, it is rather unusual since the average value of 
the nebular to pulsar emission is ( 7 ^ 7 ) ~ 4 for a large sample 
of PWNe observed with Chandra (IKargaltsev & Pavlovi!2008l) . 
The F of the nebular emission is compatible with that usually ob- 
served for PWNe (1.5 to 2.1 : lKargaltsev & Pavlovl2008l : lLi et al.l 
2008), and is also c onsistent with particle w ind models of Fermi 
shock acceleration dAchterberg et alj[2001l) . 

6.3. Morphology 

The morphology of the PWN is elongated along the east-west 
direction as reported in section l4~4l The obtained size of 5.2" of 
the X-ray nebula corresponds to a source extent of 1.5 pc at the 
distance of the SMC. The radio extent of the source is much 
larger corresponding to 21 pc with the brightest emission from 
the radio lobes, each ~ 6 pc in extent assuming the same dis¬ 
tance. This scenario is consistent with a longer cooling time for 
radio emission compared to that of X-ray emission. In addition, 
the radio nebula extends even further than the estimated SNR 
centre, consistent with the fact that the earliest electrons were 
produced when the pulsar was much closer to the SNR centre. 
Further, the long axis of the radio and X-ray nebula are aligned 
with each other. This axis of alignment points to the centre of 
the SNR, so is presumably aligned with the pulsar’s direction of 
motion. 
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6.4. Interaction of the SNR with the PWN and its stage of 
evolution 


The SNR in IKT 16 is in the adiabatic_ Sedov phase of evo¬ 
lution with an age of ~ 14.7 kyr (lOwen et al.l 1201 ll) . The 
shell of the SNR has important implications on the PWN, 
and the physics of these systems is extremely complicated due 
to th e rapid evolution of bo t h the SNR and the centra l pul¬ 
sar (ICiaensler & Slanel 120061: iGelfand et al.l 120071 l2009h . The 
evolution of a PWN inside an SNR follows three important 
evolutionary phases: an initial free-expansion in the super¬ 
nova ejecta, the collision between the PWN and SNR reverse 
shock, which crushes the PWN subjecting it to various instabil¬ 
ities, and eventually subsonic re-expansion of the PWN in the 
shock heated ejecta (iBlondin et ai ]l2001t Ivan der Swaluw et all 


2003^ Ivan der Swaluwl 120031 Van der Swaluw et ahT 2004 : 
Gelfand et al. 120071 and references therein). After the reverse 
shock interaction has died off, i.e. in the PWN re-expansion 
stage, the neutron star is often displaced from its PWN, leav¬ 
ing behind the ’’relic” PWN usually observed in radio and forms 
a ’’new” PWN comprising of freshly injected particles, and 
therefore observable in the X-rays. A later additional phase of 
”bow-shock” nebula can sometimes be identified, as the local 
sound velocity decreases with the pulsar’s progression through 
the SNR. In this phase the pulsar’s motion may become su¬ 
personic when it approaches the s hell of the SNR, and it may 
acquire a bow-shock morphology (| van der Swaluw et al.l 12003: 

Ivan der Swaluwll2003l: I Van der Swaluw et al JI2004I). 

Adopting the total mechanical energy released in the explo¬ 
sion (£51 in units of 10 51 erg) and the ambient medium density 
(no) from the adiabatic Sedov modeling in lOwen et al.1 (1201 ill , 
and assuming the mass of the ejec ta, M g , = 10 M 0 , we can c al- 
culate the reverse shock trajectory (iMcKee & Truelovell 1995l) at 
the present age of the SNR (14.7 kyr). At the distance of the 
SMC, we obtain that the reverse shock radius R s is 11 pc away 
from the centre of the SNR. This would imply that although it 
might not have encountered the smaller X-ray nebula_yet (mea¬ 
sured to be at 8 pc from the centre of the SNR iOwen et all 201 Ilk 
it would have crossed the larger radio extent of the nebula and 
interacted with the lobes. This picture is however not supported 
by the new high resolution radio data. The presence of two bright 
radio lobes more or less symmetric about the X-ray counterpart 
is consistent with an expansion of the PWN into cold ejecta in 
homologous expansion. We however point out that there may be 
uncertainties in measurement of our Sedov model parameters for 
the SNR IKT 16 (lOwen et alJl201 lh which may limit the estima¬ 
tion of R*. 


7 . Conclusions 

Thanks to the high resolution Chandra ACIS observation we 
have been able to resolve the hard X-ray emission near the cen¬ 
tre of the SNR IKT 16 into a PWN centred on a putative pulsar. 
We have imaged the PWN, constrained its geometrical parame¬ 
ters, and have measured all the spectral components separately. 
With new high resolution radio data we have determined the ra¬ 
dio morphology and its extent precisely. The main results can be 
summarized as follows: 

- The putative PWN is elongated but centred on the point 
source. Comparison with the radio counterpart indicates that 
the reverse shock may not have yet met with the PWN sur¬ 
face. 

- The point source at the centre is about three times brighter 
than the elongated feature. 


- Morphological modelling of the X-ray nebula with a PSF 
and a 2D Gaussian resulted in an FWHM of 5.2" for the 
PWN with its axis aligned with the larger radio nebula. 

- The point source at the centre has a much harder spectrum 
than the extended emission from the nebula. This points to 
the presence of a pulsar dominated by non-thermal emission. 
The expected E is ~ 10 37 erg s _I and spin period j 100 ms, 
with P ~ 10~ 13 s s -1 . This points towards a young and en¬ 
ergetic pulsar powering the PWN in IKT 16. However, the 
presence of a compact X-ray PWN unresolved by Chandra 
at the distance of the SMC cannot be ruled out. 
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